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ABSTRACT 

We report the first constraints on the properties of weakly interacting low-mass dark matter (DM) 
particles using asteroseismology. The additional energy transport mechanism due to accumulated 
asymmetric DM particles modifies the central temperature and density of low-mass stars and sup- 
presses the convective core expected in 1.1-1.3 M Q stars even for an environmental DM density as low 
as the expected in the solar neighborhood. An asteroseismic modeling of the stars KIC 8006161, HD 
52265 and a Cen B revealed small frequency separations significantly deviated from the observations, 
leading to the exclusion of a region of the DM parameter space mass versus spin-dependent DM-proton 
scattering cross section comparable with present experimental constraints. 

Subject headings: asteroseismology - dark matter - stars: individual (a Cen B, KIC 8006161, HD 
52265) 



1. INTRODUCTION 

The identification of the nature of the dark matter 
(DM) of the Universe i s a major open problem in modern 
Among the diverse strategies for 
y of the possible impact of DM in 



physics 
DM sear* 



( |Bertone|2010| . 

rches, the study 



as a complementary approach to the DM problem 


Spol- 


yar et al. 2008| |Scott et al 


2009 Casanellas & 


Lopes 


2009| IZackrissonet al.|2010; 


Sivertsson & Gondolo 


2011 


Casanellas & Lopes 2011a 


Scott et al.||2011 Li et al. 


2012 Hie et ai.||2012| |C6rsico et al.||2012). In partic- 



sic matter-antimatter asymmetry dKaplan et al. 


2009 


Davoudiasl et al.|2011 


Blennow et al.|2012 


) do not anni- 



hilate after gravitational capture by compact astrophys- 
ical objects and can therefore strongly influence their in- 



ternal structure ( Griest fe Seckel||1987[ ). Thus, both the 



observation or the lack ot observation of the impact of 
asymmetric DM (ADM) on the properties of stars can 
be used to put constraints on the characteristics of these 
DM candidates. 

The interior of the Sun, being known with a high ac- 
curacy thanks to solar neutrinos and helioseismic data, 
is an excellent laboratory to probe the existence and the 
properties of ADM particles. Such particles remove en- 
ergy from the inner ~4% of the Sun, leading to a reduc- 



isothermal core (|Taoso et al. 


2010 


Frandsen & Sarkar 


2010 Lopes & Silk 2012). In particular, ADM candidates 



scattering cross sections may influence the internal solar 
structure so strongly that they would produce clear sig- 
natures in the low-de gree frequency spacings and in the 
solar gravity modes (iLopes fc Silk 2010; |Cumberbatch 



et al.|2010| |Turck-Chieze et al.|2012f ~Tnterestingly, low- 



mass weakly interacting massive particles (WIMPs) with 
similar characteristics provide an explanation for the sig- 
nals in various direct detection experiments, strengthen- 
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ing the motivation for the search of indirect signatures 
of these particles. 

It has also been shown that these low-mass ADM 
candidates may produce mar ked effects in very low - 
mass stars and brown dwarfs ( Zentner fe Hearin||2011 ). 
In environments with high ADM densities, solar-like 
stars may sho w significant devia tions in their evolu- 
tionary tracks (Iocco et al. 2012). Also neutron stars, 



due to their compactness, capture DM very efficiently 
and m ay be strongly influenced by the accumulation of 



ADM dBertone fe Fairb airn 
120111 iLeung et al.||2012 



2008 Kouvaris & Tinyakov 



ere we will show that, even 
for a DM density as low as the e xpected in the sola r 
neighborhood, p x = 0.4 GeV cm 



( Garbari et al.|2012 ) 

main-sequence stars with masses similar to that ol the 
Sun present distinct signatures of the captured ADM. 

With the advent of asteroseismology, a precious in- 
sight into the stellar interi ors is nowadays po ssible for 
the fi rst time. The CoR oT ( |Michel et ail2008| ) and Ke- 
pler ( Bruntt et al.||2012| missi ons have already de tected 
oscillations m about 500 stars (Chaplin et al. 2011 1. This 
fact has allowed to test theories of stellar evolution and 
to pr obe the stellar cores with an unprece dented preci- 
sion ( Garcia et al.|2010 Bonaca et al. 2012). The seismic 
analysis of stars other than the Sun is complementary to 
helioseismic DM searches because it allows the study of 
stars with lower masses, which are more strongly influ- 
enced by DM, and stars whose dominant energy trans- 
port mechanisms may change due to the DM influence. 
In this Letter we will demonstrate, by studying the case 
of the stars KIC 8006161, HD 52265 and a Cen B, that 
present asteroseismic observations do constrain a signif- 
icant region of the DM parameter space. 

2. INTERACTION OF DARK MATTER AND STARS 

Nearby stars are embedded within the halo of DM par- 
ticles that is presently believed to permeate our Galaxy. 
If these DM particles have a non-negligible scattering 
cross section off baryons (so they are WIMPs), then 
some of them may collide with the nucleons of the stel- 
lar plasma, losing kinetic energy. A fraction of these 
DM particles is gravitationally captured by the stars. To 
calculate the capture rate we follow the formalism that 
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TABLE 1 

Constraints on the stellar characteristics adopted for the modeling and selected results. 



Star 


M (M ) 


R (R Q ) 


L (L Q ) 


Teff (K) 


(Z/X) 3 


(Ai/ n , ) a (a*Hz) 


<5i/ 2> a (AtHz) 


IVH • OUuUlUl 
















Observations' 5 


0.92-1.10 


0.90-0.97 


0.61 ± 0.02 


5340 ± 70 


0.043 ± 0.007 


148.94 ±0.13 


10.10 ±0.16 


Stand, modeling 


0.92 


0.92 


0.63 


5379 


0.039 


149.03 


10.12 


DM modeling 


0.92 


0.92 


0.63 


5379 


0.039 


149.08 


9.13 


HD 52265 
















Observations' 5 


1.18-1.25 


1.19-1.30 


2.09 ±0.24 


6100 ± 60 


0.028 ± 0.003 


98.07 ±0.19 


8.18 ±0.28 


Stand, modeling 


1.18 


1.30 


2.22 


6170 


0.028 


97.92 


8.16 


DM modeling 


1.18 


1.30 


2.22 


6170 


0.028 


98.05 


7.65 


a Cen B 
















Observations' 5 


0.934 ± 0.006 


0.863 ± 0.005 


0.50 ±0.02 


5260 ± 50 


0.032 ±0.002 


161.85 ± 0.74 


10.94 ± 0.84 


Stand, modeling 


0.934 


0.868 


0.51 


5245 


0.031 


162.56 


10.23 


DM modeling 


0.934 


0.868 


0.51 


5230 


0.031 


162.45 


8.95 



a Averages for the intervals 2750 < v(iiRz)< 3900 (KIC 8006161), 1600 < i/(/jHz)< 2600 (HD 52265), and 3300 < v(hRt,)< 5500 
(a Cen B). 
b Data, from 



Mathur et al 
|et al.||2005] (lor a Cen B 



(2012 i and 



Bruntt et al. 



120121 (for KIC 8006161), 



Ballot et al 



(20111 (for HD 52265), and Kjeldsen 



m x = 5 GeV, <t x> sd = 3 x W~ ib cm 2 , p x = 0.4 GeV cm" 



was early deve l oped b y Gould (1987), as implemented in 
Gondolo et al. ( 2004 1 . We assume a Maxwell-Boltzmann 
distribution of the velocities of the DM particles, with 
a dispersion v x — 270 km s^ 1 , and a stellar velocity of 
v+ = 220 km s _1 . The expected deviation from the men- 
tioned fiducial values for the specific stars studied in this 
work may lead to a maxi mum error on the ca pture rate of 
approximately 15% (see Lopes et al. (2011) for details). 

In the asymmetric WIMP scenario, the annihilation 
cross sections required to match the DM relic density are 
larger than in the standard WIMP scenario, dep ending 
on the degree of asymmetry (Iminniyaz et al. 2011 1. Nev- 
ertheless, in the asymmetric scenario the particles con- 
centrated on the stellar core cannot find a partner with 
whom to annihilate (the DM particle is not its own anti- 
particle) and therefore their number grows indefinitely 
while more particles are being captured. The same effi- 
ciency in the stellar accumulation of DM occurs for very 
feebly annihilating ((cr a v) < 10 -33 cm 3 s _1 ) Majorana 
DM particles. In other cases, the captured DM would 
have no relevant impact on nearby low-mass stars. 

The evaporation of DM particles can be neglected for 
the ADM candidates and stars considered in this work: 
while a low stellar mass tends to favor evaporation, this 
fact is compensated b y the cooler stellar tem peratures, 
Evap oc e -GMm x /RT flGriest fc Seckel 
ment with the results of |Zentner & Mearm 
~ 0.1 M stars. 

The DM particles captured in the stellar core provide 
a new energy transport mechanism that removes energy 
from the center of the star. The efficiency of this mech- 
anism depends mainly on the ratio between the mean 
free path of the WIMPs inside the stellar plasma l x and 
the characteristic radius of the WIMPs distribution in 
the core of the star r v 



1987|) , in ag ree- 
pJllT for 



(Gilliland et al. 



1986). For most 



of the WIMP-proton SD scattering cross sections ct x ,sd 
considered here, l x > r x and the energy transport by 
WIMPs is non-local. On the other hand, for large values 
°f a x,SD> < r x so t ne WIMPs are in local thermal equi- 
librium with the baryons. The latter regime applies only 
to values of <j x .sd which are not considered in this work 
(<?x,S D ~ 10 -33 cm 2 ). However, we follow the prescrip- 



tion described in Gould fe Raffelt| ( 1990 ) that extends the 
formalism developed for the local thermal equilibrium to 
other regimes by the use of tabulated suppression factors. 

The DM capture and e nergy transp ort mechanisms 
were implemented in CESAM (|Morel 1997), a sophisticated 
stellar evolution code. In the case of the Sun, the result s 



of our modified solar model (e.g. Lopes & Silk (2012)) 



are i n agreement with those of other codes in the litera- 



ture ( |Taoso et al"]|2010| |Frandsen fc Sarkar||2010| ). The 

observational constraints used tor the modeling of the 
stars KIC 8006161, HD 52265, and a Cen B, as well 
as the results of some selected models with and without 
taking into account the DM effects, are summarized in 
Table [TJ 

3. IMPACT OF ADM ON THE PROPERTIES OF LOW-MASS 

STARS 

3.1. Modifications of central temperature and density 

The main signature of the additional DM cooling mech- 
anism is a decrease in the central temperature and an in- 
crease in the central density. These variations are shown 
in Figure [l] for several DM-modified stellar models, cal- 
ibrated to reproduce the observed properties of the star 
KIC 8006161, for a range of DM masses and SD scatter- 
ing cross sections. Compared with the standard model- 
ing, for m x = 5 GeV and c x ,sd = 3x 10~ 36 cm 2 we 
found a ~ 9% decline in the central temperature. The 
variations on the internal properti es are larger than t hose 
reported in the case of the Sun (Taoso et al. 2010) be- 
cause the importance of the energy transported by the 
WIMPs (e x . trans oc C x oc over the thermonuclear 
energy (e nuc i oc M 3 5 ) increases when the stellar mass 
decreases. In particular, in our computations we found 
the DM cooling to reduce the T c of 0.7 M stars nine 
times more efficiently than for 1.1 M stars. This fact 
reinforces the potential advantages of performing DM 
searches in stars other than the Sun. 

3.2. Suppression of convective core 

In the standard picture of stellar evolution, stars with 
masses greater than 1.1 M are expected to keep a con- 
vective core during most of the main sequence, while 
stars with lower masses quickly lose their convective 
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Fig. 1. — Central temperatures (top) and densities (bottom) of 
the DM-modified stellar models that reproduce the observed prop- 
erties of the star KIC 8006161. 

cores. Convection arises when the gradient of temper- 
ature in the core is so steep that a rising bubble of 
plasma does not cool enough with its adiabatic expan- 
sion, so that it continues to rise, leading to a convective 
instability. If the temperature gradient is reduced by 
an additional mechanism such as the energy transport 
by WIMPs, then the conditions for convection may no 
lo nger be achieve d. This possibility was first suggested 




Renzinij (1987), where the suppression of convection 
in horizontal branch stars was predicted using analyti- 
cal approximations. This scenario must not be confused 
with the creation of an unexpected convective core in 1 
M Q stars due to the self-annihilation of DM pa rticles cap- 
tured in halos w ith very high DM densities ( |Casanellas| 
fe Lopes||2011b| ). 

The reduction of the temperature gradient in the stel- 
lar interior due to the additional cooling by WIMPs was 
found to suppress the convective core expected in stars 
with masses slightly greater than that of the Sun. The 
standard modeling of the star HD 52265 predicted a con- 
vective core during all the main sequence, but this con- 
vective core rapidly disappeared when the energy trans- 
port by WIMPs was taken into account (see Figure^ a)). 
The range of DM masses and SD scattering cross sections 
for which the suppression of the HD 52265 convective 
core is expected is shown in Figure |2^b)). Interestingly, 
hints of the sign atures of a convectiv e core in HD 52265 



were reported in Ballot et al 
elusive information can 



However, no con- 
e extracted until there is no 



definit ive diagnostic of it s presence or its absence (see 
also|Escobar et aL (2012)). 



4. ASTEROSEISMIC DIAGNOSTIC OF THE PRESENCE OF 

DM 

The characteristic signatures reported in the last sec- 
tion are potentially detectable with the analysis of the 
stellar oscillations. Asteroscismology is presently show- 
ing its power in determining with high precision not 
only the global properties of stars but also their inter- 
nal structure. In particular, the small frequency sepa- 
rations of low angular degree (I = 0) and radial order 
n, Si>q2 = v n .o — v n ~i.2, have been shown t o provide use - 
ful information about the core of the stars ( |Gough|1986[ ). 
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Fig. 2. — (a) Size and duration of the convective core in the 
modeling of the star HD 52265 in the classical picture (gray) and 
taking into account the energy transport due to ADM particles 
with m x = 5 GeV and <r x , SD = 1.5 X 10" 36 cm 2 (blue), (b) The 
presence of a convective core in HD 52265 depends on the mass 
and SD scattering cross section of the DM particles. 

Thus, we would expect the seismic parameter (61/02) to be 
sensitive to the modifications introduced by the WIMPs 
on stars. 

We have computed the oscillation frequencies and 
separations of the DM-modified stellar models of KIC 
8006 161, HD 52265 and a Cen B using the ADIPLS pack- 
age ( Christensen-Dalsgaard 2008). In order to disentan- 
gle the effects ot DM trom those arising from the variation 
of the stellar parameters, a very precise determination 
of the latter is of utmost importance. Although astero- 
scismology has already provided very accurate measure- 
ments of t he mass and ra dius of KIC 8006161 (Mathur 
let al.|2012|) and HD 52265 " 



( |Escobar et al.|2012 ), with un- 
certainties of the order of 1%, we preferred to focus here 
on the case of a Cen B, a star whose fundamental pa- 
rameters are independently measured with high precision 
(see Table T| thanks to its proximity and because it be- 
longs to a binary system. The mass of this star has been 
determined from the radial velocities of a Cen A and B, 
its effective temperature from high quality spectra, its lu- 
minosity from photometric data and its radius from mea- 
sureme nts of its angular diam eter combined with paral- 
lax (see Kjeldsen et al. ( 2005 ) and references therein) . 

All stellar models used to create Figure [3] reproduce 
the measured M, L, R, T e rf, (Z/X) s and mean large 
frequency separation (Av n i) of a Cen B within the ob- 
servational error. However, while models without DM 
are also able to reproduce the observed mean small fre- 
quency separation (^^02), we found that the stellar mod- 
els with a strong influence of DM predict a (Sv^) sig- 
nificantly deviated from the observed value. The black 
lines in Figure [3j labeled 2a and 5er, show the DM char- 
acteristics corresponding to the calibrated models that 
predicted a ((^02) with a difference of two and five times 
the observational error, respectively, from the observed 
value. The dashed black lines around the 2a line show 
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m x (GeV) 

Fig. 3. — Deviation of the small frequency separation (5uq2) of the 
DM-modificd stellar models from the true value measured in a Cen 
B. All the stellar models are calibrated to fit the M, R, L, T e ff, 
(Z/X) s and {Av n j) of a Cen B within the observational error. The 
dashed black lines around the 2cr line show the uncertainty in (<Wo2) 
arising from the observational error in the stellar characteristics. 



the uncertainty in the modeling when the observational 
errors in the stellar characteristics M, L, R, T e ff and 
(Z/X) s are taken into account. This uncertainty cor- 
responds to the standard deviation on ((^02), evaluated 
computing 2600 valid models of a Cen B and including 
also the uncertainty in the capture rate C x from varia- 
tions in the DM halo parameters and the stellar velocity, 
as discussed in Section [2] The dashed lines around the 
5cr line are not shown for clarity, but they would appear 
narrower because {8v 02 ) varies more abruptly in that re- 
gion of the plot. Therefore, we conclude that present 
asteroseismic measurements of a Cen B disfavor the ex- 
istence of DM particles with parameters above the 2a 
line with 95% confidence level. 

Similarly, the presence of a convective core leads to 
strong asteroseismic signatures. The mixing of elements 
in convective regions introduces sharp structural varia- 
tions in the border with radiative regions that produce 
a clear oscillatory signal in the frequency spectrum. It 
has been shown that this feature may be used to detect 
and measure the size of a convective core thro ugh aster- 
oseismic parameters such as r m , r ln or dr 2i3 (|Cunha & 



; ps 

Brandaop011| |Silva Aguirre et al.||2011|) . If these aster 



oseismic diagnostic tools succeed in the confirmation of 
the presence or the absence of a convective core in a star 
with 1.1-1.3 Mq, this hypothetical measurement may be 
used to place further constraints on the nature of the DM 
particles. The characteristic and localised effects of DM 
should allow the disentanglement of its signatures from 
standard processes. Remarkably, several stars with the 
appropriate characteristics are presently being observed 
by the C0R0T and Kepler missions. 



5. CONCLUSIONS AND DISCUSSION 

We have shown the strong signatures that asymmetric 
DM particles with low masses and large SD scattering 
cross sections with baryons produce on low-mass stars. 
We have focused in the study of the stars KIC 8006161, 
HD 52265 and a Cen B, revealing large modifications 
in the central temperatures and densities of the models 
and the suppression of the convective core expected in 
1.1-1.3 M stars. 

In the case of a Cen B, we have shown that the 
asteroseismic parameter (SV02) can be used to impose 
competitive constraints to the DM characteristics. In 
particular, we were able to exclude with 95% confi- 
dence level ADM candidates with m x ~ 5 GeV and 
a x,SD > 3 x 10~ 36 cm 2 . These new constraints are com- 
parable with the present limits from direct detection ex- 



per iments (oy.sp 2 x 10 



for m, 



5 GeV, 



— ^ ~ Y — — - 1 * ' 

Archambault et al. (20121)) because the sensitivity of 



see 

the detectors drops at low WIMP masses 

Interestingly, low-mass WIMPs with similar character- 
istics have been advocated to explain the signals in the 
DAM A/LIBRA and CoGeNT detectors in terms of SD 
collisions. In ADM models the low mass of the WIMPs 
is strongly motivated because the relic density of DM is 
determined by the baryon asymm etry of the Universe , 
leading to Qdm ~ (niDM /mb)Slb (Kaplan et al. 2009). 
Our approach may provide a complementary test ot these 
low-mass WIMP models, in the context of controversy 
over the incompatible results between different direct de- 
tection experiments. 

Asteroseismology thus arises as a promising strategy 
for indirect DM searches. Compared to helioseismology, 
the asteroseismic searches of DM allow the study of stars 
with masses lower than that of the Sun, which are more 
strongly influenced by the additional cooling mechanism 
provided by the DM particles. In addition, the aster- 
oseismic confirmation of the presence or the absence of 
convective cores in 1.1-1.3 M Q stars, such as HD 52265, 
may provide further constraints on the nature of DM. 

The future perspectives of this approach are also excit- 
ing. If the small frequency spacings are identified in the 
oscillations of stars located in environments with high 
expected DM densities, such as globular clusters, then 
the sensitivity of the approach proposed in this work 
will reach much smaller WIMP-proton scattering cross 
sections and larger WIMP masses. Moreover, in the 
event of a successful identification of the properties of 
DM after hypothetical positive results in different exper- 
iments, asteroseismology may allow the determination of 
the density of DM at any specific location where a star 
is observed. 
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